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Abstract. We consider the Schrodinger operator —A + q in domains of the form 
R = {x € W 1 : < Xi < a,i, i = l,..,n} with either Dirichlet or Neumann 
boundary conditions on the faces of R, and study the constraints on q imposed by 
fixing the spectrum of —A + q with these boundary conditions. We work in the 
space of potentials, q, which become real-analytic on M n when they are extended 
evenly across the coordinate planes and then periodically. Our results have the 
corollary that there are no continuous isospectral deformations for these operators 
within that class of potentials. This work is based on new formulas for the trace 
of the wave group in this setting. In addition to the inverse spectral results these 
formulas lead to asymptotic expansions for the traces of the wave and heat kernels 
on rectangular domains. 

§1. Introduction 

We consider the Schrodinger operator, H = — A+q(x), in the rectangular domain 
R = {x e R n : < Xi < di, i = 1, .., n} and on each of the 2™ faces of dR, Rn{xi = 
or Xi = a>i} we impose either Dirichlet or Neumann boundary conditions. Let 

Ml < fJ-2 < fJ-3 < ■ ■ ■ 

be the spectrum of the self-adjoint operator, also denoted by H, that is obtained in 
this way. We take R and as given, and study what constraints these data 

put on q. 

This work depends strongly on the geometry of rectangular domains, and involves 
many reflections, translations, etc. Since all of this will be easier to follow once one 
has seen it in a simple case, we have moved the statements and proofs of the general 
results to §5. In this Introduction and §2-§4 we will restrict ourselves to the case 
n = 2 with Dirichlet boundary conditions, i.e. we will work in the rectangle 

R = {(xi, X2) '■ < x\ < a, < X2 < b}, 

with the boundary condition u = on dR. 

For our results we need to extend q to an even periodic potential on M 2 : we 
extend q to Q on M 2 by defining 

Q(-x 1 ,x 2 ) = q(xi,x 2 ) for x € R, 

Q(x\, —x 2 ) = Q(x\, x 2 ) for \x\\ < a, < x 2 < b, and 

Q{x\ + 2ma, x 2 + 2nb) = Q(xi,x 2 ) for < a, \x 2 \ < b, and (m, n) € Z 2 . 

Thus Q is periodic with respect to the lattice L = {(2ma,2nb), m, n G Z}. To L 
we associate the dual lattice 



L* = {6 G R 2 : 5 ■ d € Z for all d € L}, 
1 
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and expand Q in a Fourier series 

Q(x) = a & e 2 ^\ 
<5eL* 

Let S denote the set of elements of S which are maximal in the sense that {5 ■ d, d G 
L} = Z. For our lattice S = {(m/(2a),n/(2b)) : m and n are relatively prime }. 
Then we can decompose Q 

oo 

5eS k=-oo 8€S 

where 

oo 

Qs(s)= a kse 2mks - 

k = — oo 

Note that this decomposition holds if, but only if, a( 0j o) = J R qdx = 0. How- 
ever, in all dimensions (see §6) one can recover j R qdx from the asymptotics of 
the trace ^]exp(— /j,jt) as t — > 0. So given isospectral potentials, subtracting this 
constant from each of them, one can replace them by isospectral potentials satis- 
fying J R qdx = 0. Note also that Q_s(s) = Qs{— s) which gives rise to the factor 
1/2 in these formulas. As in [ERT1] we call Qs a "directional potential". Since 
Q(x) = Q(—x), we have Qs(s) = Qs(—s). The directional potentials corresponding 
to coordinate directions, i.e. <5i = (l/(2a),0) and 62 = (0, 1/(26)), only depend on 
the coordinates x\ and X2 respectively, and we denote them by q\{x\) and (72(^2)- 
The main result of the first part of this article is: 

Theorem 1.1. Assume that —A + q and — A + g have the same Dirichlet spectrum 
on R. If a 2 /b 2 is irrational and the extensions Q and Q of q and q described above 
are real-analytic on M. 2 , then 

a) for any 5 6 S with no zero components the operators — \5\ 2 d 2 /ds 2 + Qs(s) and 
— \S\ 2 d 2 /ds 2 + Qs(s) with periodic boundary conditions on [0,1] have the same 
spectrum. 

b) the pairs of operators, —d 2 /dx 2 + qi(x\) and —d 2 /dx\ +qi(x\), and, —d 2 /dx 2 + 
92(^2) and —d 2 /dx 2 + (72(^2), have the same spectrum for Dirichlet boundary con- 
ditions on [0, a] and [0.6] respectively. 

Part a) of Theorem 1.1 was a surprise. In the case of periodic boundary condi- 
tions it was natural that the directional potentials would be isospectral for periodic 
boundary conditions, but we did not anticipate that this would be true for Dirich- 
let conditions. Another unexpected result in comparison to the periodic case is 
the rigidity of the directional potentials Qs for 5 in the "oblique directions" in a) . 
Since Qs is even, the isospectral set for — \5\ 2 d 2 / ds 2 + Q s(s) with periodic boundary 
conditions on [0, 1] is either finite or a Cantor set, depending on whether a finite or 
an infinite number of "gaps" in the spectrum are open (see [GT]). 

While part b) is the most one can say for pairs of isospectral potentials of the 
form q(x) = qi(xi) + 92(^2)) i- e - f° r potentials of that form b) implies that q 
and q are isospectral on R, this is not the end of the story. Since for our class of 
potentials, the reduced potentials qi(x\) and (72(2:2) also satisfy the one-dimensional 
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version of the symmetry conditions, it turns their Dirichlet spectra on [0, a] and [0, b] 
determine their periodic spectra on [—a, a] and [—b, b] respectively. Thus, like the 
directional potentials in the oblique directions, the set of isospectral potentials in 
the coordinate directions is discrete. In other words within our class of potentials q 
the set of isospectral potentials for the Dirichlet condition on a rectangle is discrete. 

Theorem 1.2. If o? jb 2 is irrational, then any continuous curve q t of isospectral 
potentials such that the extensions Q t are real analytic for all t must be constant, 
q t = qo for all t. 

The proof of Theorem 1.1 is given in §2 - §4. It is an extension of the proof 
of the analogous result, [ERT1, Theorems 6.1 and 6.2], for periodic boundary con- 
ditions. The new ingredients in the proof are in the analysis of the trace of the 
fundamental solution of the wave equation in R with the potential q. The Dirichlet 
trace is significantly different from the trace for periodic boundary conditions, but 
the new terms in the Dirichlet trace "telescope" in a way which simplifies their 
contribution to the singularities of the trace. The resulting formulas reveal close 
relations between the Dirichlet and periodic traces, and they also make it possi- 
ble to compute the singularities in the trace in a way that identifies contributions 
with the underlying geometry. For this we found it simpler to use the classical 
Hadamard construction of parametrices for the wave equation instead of Fourier 
integral operators. Given the lemma that isospectral potentials in the coordinate 
directions form a discrete set, Theorem 1.2 is an immediate corollary of Theorem 
1.1. The lemma is proven in §5 (Lemma 5.3). 

In [ERT2] we showed that, under the assumptions in Theorem 1.1, further anal- 
ysis of the heat trace associated with this problem reduced the set of isospectral 
potentials for —A + q with periodic boundary conditions to just q(x) and q(—x) for 
many choices of q. We plan to investigate the applicability of that analysis to the 
Dirichlet problem in future work. 

The following result does not require the analyticity of the extended potential 
Q or the irrationality of a 2 jb 2 : 

Theorem 1.3. Suppose that q and q have the same Dirichlet spectrum on R, and 
have extensions Q and Q in C°°(IR 2 ). Then, if q(x) = qi(xi) + q2 (#2), there are 
smooth potentials q\ and q\ such that q(x) = qi(x\) + 52(^2)- 

The proof of Theorem 1.3 follows the outline of the proof of the corresponding 
result, Theorem 4.1, in [ERT1]. The authors and E. Trubowitz proved a version of 
Theorem 1.3 in 1981, using the approximate eigenfunction construction in [ERT1, 
Section 3(b)], and we still plan to return to approximate eigenfunctions in future 
work. 

As we said earlier, the generalization of Theorems 1.1 and 1.2 to rectangular 
domains in M n with arbitrary assignment of Dirichlet and Neumann boundary con- 
ditions on their faces is given in §5. The constructions and arguments are sufficiently 
analogous to the two dimensional Dirichlet case that they can be presently rather 
concisely. Section 6 is devoted to the expansion of the singularities of the wave trace 
at t = |<Z|, d € L\0. This is used for the proof of Theorem 1.3 and its generalization 
to higher dimensions, Theorem 6.1. 

In §7 we use the Hadamard expansion of the fundamental solution to expand 
the wave trace near t = in distributions of increasing orders of homogeneity. This 
expansion can be used to derive asymptotics for the heat trace, ^]exp(— Hjt), as 
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t — > 0. While there is a huge literature on the asymptotics of heat traces, the 
relatively simple special case we treat here appears to be new. The expansions 
reflect the singularities in the rectangular geometry, and contain terms which are 
absent for smooth boundaries. For instance, in dimension two with the Dirichlet 
condition f R q 2 dx no longer appears as a spectral invariant, being replaced by 
J R q 2 dx — (tt/2) q(Pi), where Pi, i = 1, ..,4, are the corners of the rectangle. 

There is an extensive literature on multi-dimensional inverse spectral problems, 
and we have included a sampling in the references. In the setting of inverse spectral 
problems for —A + q on given domains, the literature includes Guillemin-Kazdan 
[GuK], Guillemin [Gul] and [Gu2], Gordon-Kappeler [GK] and, recently, Gordon- 
Schiith [GS]. 

§2. Fundamental Solutions, Traces and Cancellations 

In this section it suffices to assume that q extends to a sufficiently smooth Q £ R 2 . 
Let E(t, x, y) be the fundamental solution for the initial value problem 

u tt = Au-Qu in R t x R 2 , u(0, x) = f(x), u t (0, x) = 0. (1) 

It follows from the properties of the wave front sets of solutions of the wave equation 
(see, for example, [H61]) that E(t,x,y) is an even distribution in t on R, depending 
smoothly on x and y. Here and elsewhere we follow the convention of writing 
distributions as functions, because we think it makes the manipulations easy to 
follow. 

Using the lattice L from Section 1, the fundamental solution for (1) in R with 
Dirichlet boundary conditions on dR can be written 

D(t,x,y) = 

y~][E(t,d + x, y) - E(t, d 1 -x 1 ,d 2 +x 2 ,y)-E(t, di + x 1 ,d 2 -x 2 ,y) + E(t,d-x,y)]. 

(2) 

This construction is just the "Method of Reflection". To check it, note first that 
for (x,y) £ R x R the sum is finite by finite speed of propagation for the wave 
equation. If x\ = or x 2 = each term in the sum vanishes. If x\ = a, then 
±E(t,2ma + a,2nb ± x 2 ,y) cancels ^E(t,2(m + l)a — a,2nb ± x 2 ,y), and the 
analogous cancelations occur for x 2 = b. Since E(0,x,y) = 5(x — y), all terms 
in the sum with d ^ vanish for (x,y) in the interior of R x R, and we have 
D(0,x,y) = 5(x — y) and dtD(0,x,y) = 0. The symmetries of Q imply that all 
terms E(t, d\ ± x±, d 2 ± x 2 ,y) are solutions of Uu — A x -u + qu = for x £ R. 

Note that D(t, x, y) like E(t, x, y) is a distribution in t on R depending smoothly 
on x and y, since each of the terms in (2) has this property, and only finitely 
many of these terms are nonzero when t lies in a bounded interval. Therefore the 
trace Tr(D{t)) = J R D(t,x,x)dx exists as a distribution in t on R. Note that 
f_ oo D(t,x,y)p(t)dt,p £ Cq°, is the kernel of a trace class operator on L 2 (R). 
Therefore we have the identity 

IY][p(VJ^) + P(-VJ^)} = / (/ D(t,x,x)p(t)dt)dx 

z JR JR 
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or more concisely in the sense of distributions on R 

Vcos( v ^t)= / D(t,x,x)dx = Tr D(t). (3) 
Jr 

To exploit some symmetries it is convenient to use the rectangle 

-Ro = {(xi,x 2 ) ■■ \xi\ < a, \x 2 \ < b}. 
The symmetries of E(t, x, y) corresponding to the symmetries of Q(x) imply that 

4Tr(D(t))= D(t,x,x)dx. (4) 

J R n 



^From (3) and (4) one sees that the spectrum of H on R determines J Rq D(t, x, x)dx.| 
Expanding this using (2), we have 

/ D(t,x,x)dx = y~^[ / E(t,d + x,x)dx — / E(t, d\ — xi,d 2 + x 2 ,x)dx— 

J Ro d^L 

/ E(t,d 1 +x 1 ,d 2 -x 2 ,x)dx+ / E(t,d-x,x)dx] = def e ++ (t)-e~ + (t)-e + ~ (t)+e (t)M 
JRo JRo 1 

There are cancellations in e~ + (t), (t) and e~~(t) which reduce them to sums 
of fewer terms. We will derive this reduction for e (t); the reductions for e + ~(t) 
and e '"(t) are essentially corollaries of that case. 

Since the periodicity of Q implies E(t, d + x,d + y) = E(t, x, y) for all d € L, 
substituting x\ + 2a for x\ gives 

/■OO 

/ E(t, d — x, x)dx\ = 

J a 

/oc poo 
E(t,d-(2a,0)-x,(2a,0)+x)dx 1 = E(t, d - (4a, 0) - x, x)dx x . (5) 
-a J —a 



Likewise 



/ — a pa 
E(t,d- x,x)dxi = E(t,d+ (4a, 0) - x,x)dxi. (6) 
-oo J —oo 



Since E(t, x,y) =0 when \x — y\ > \t\ 



/O nOO 
/ E(t, (-4ma,d 2 ) - x,x)dx = 0, (7) 
-b J -a 

when |(— 4ma, d 2 ) — 2x\ > \t\ for \x 2 \ < b and x\ > —a. Thus, for any R there is an 
M(R) such that (7) holds for m > M(R) when \t\ < R. Thus, writing 

/a poo poo 

E(t,d — x,x)dx\ = / E(t,d — x,x)dx 1 — / E(t, d — x, x)dx\, 
-a J —a J a 



(5) yields 

/a fco foo 

E(t,d — x,x)dxi = / E(t,d — x,x)dx 1 — / E(t, d — (4a, 0) — x, x)dxi, 
-a J —a J —a 



and we have the following telescoping sum formula 

oo r b r-a 



^ fO pa 

/ / E(t,(—2ma,d 2 ) — x,x)dx = 

m=l 

/b fOO fb poo 

/ E(t, (—2a, d 2 ) — x,x)dx + / / E(t,(—4a,d 2 ) — x,x)dx 
-b J — a J — b J — a 



/b poo pb f'OO 

I E(t,(a,d 2 )-x,(a,0) + x)dx + / E(t, (0,d 2 ) - x,x)dx, (8) 
-bJO J-bJa 

where to get the last equality we replace x\ by —a+x\ in the first integral and x\ by 
— 2a+x\ in the second integral and use the periodicity E(t, (2a, 0) — x, (2a, 0) + x) = 
E(t, —x,x). 

If one begins with 

/a pa f — a 

E(t,d — x,x)dx\ = / E(t,d — x,x)dx\ — I E(t, d — x, x)dx\, 
-a J — oo J — oo 

and uses (6), the same reasoning leads to 

°° fb na 

/ / E(t, (2ma,d 2 ) — x,x)dx = 

m=1 J-bJ-a 

/b pa pb pa 

I E(t, (2a, d 2 ) - x, x)dx + / / E(t, {^a, d 2 ) - x, x)dx 
-b J — oo J —b J — oo 

/b pO pb p — a 

/ E(t,(a,d 2 )-x,(a,0)+x)dx+ / E(t, (0,d 2 ) - x,x)dx, (9) 
-hi/— oo J — b J — oo 



where to get the last equality we change x\ to a + xi in the first integral, change 
X\ to xi + 2o in the second integral and use the periodicity in the second integral. 
Combining (8) and (9), we have 



00 fb fa 

/ / E(t, (2ma, d 2 ) — x, x)dx 

= _„ J -b J -a 



m= — oo 



/b f 00 />fo /-oo 

/ E(t,(0,d 2 ) - x,x)dx+ / (a,d 2 ) - x, (a,0) + x)dx (10) 

-b J — 00 J —b J —00 

Now we can apply the same reasoning in the second variable. That gives 
e (*)=defXl/ E(t,d-x,x)dx = / E(t,-x,x)dx+ E(t,(0,b)-x,(0,b)+x)da 



+ / E(t, (a,0) -x,(a,0) +x)dx+ I E(t,(a,b) - x,(a,b) + x)dx. (11) 
jr? Jr 2 

For e '"(t) and e" 1 (i) there are only cancellations in the sums in the first and 
second components of the lattice vectors respectively, resulting in the formulas 

e '"(t) = [/ (/ E(t,—xi,2nb + x 2 ,x)dx 2 )dxi 

/oo /-b 
(/ a- x 1 ,2nb + X2,a + x 1 ,X2)dx 2 )dx 1 ], (12) 

-00 J — b 

and 

e H (i) = [/ (/ E(t,2na + x\, —X2, x)dx 1 )dx2 

/oo />a 
(/ J5(t, 2na + x 1 , b - x 2 ,xi, b + x 2 )dxi)dx 2 ]. (13) 
-00 J —a 

Since E(t,x,y) = when |x — y| > t the integrals over M 1 and M 2 in (11), (12), 
(13) can be written with finite limits of integration. For example 

/b poo 
I E(t, a — x, 2nb + x 2 , a + X\, x 2 )dxidx 2 
-b J — 00 

= / / E(t, a — Xi, 2nb + x 2 , a + x\, X2)dxidx2 and 

J-b J4x 2 <t 2 -4,n 2 b 2 

/ E(t, (a, b) — x, (a, b) + x)dx = / E(t, (a, b) — x, (a, b) + x)dx. 

Jr? J4\x\ 2 <t 2 

Analogous formulas hold for other terms in (11), (12), (13). Moreover, for each 
term in Tr(D(t)) there is a t c such that domain of dependence considerations, i.e. 
E(t, x, y) = for \t\ > \x — y\, imply that the term vanishes for \t\ < t c . In the next 
section we will see that these terms are real-analytic in t when \t\ > t c . 

We can rewrite Tr(D(t)) in a form which has a geometric interpretation. We 
will write 

TrD(t)= I D(t,X,x)dx = J2 D rnn(t), 

where D mn (t) will correspond to the lattice point (2ma,2nb) £ L and have 

t c = 2\J m 2 a 2 + n 2 b 2 . When mn 7^ 0, D mn (t) is simply the following term from 

e++ 



D mn (t) = / E(t, x + (2ma, 2nb), x)dx, 
JRo 

Note that these terms appear in the wave trace in the case of periodic boundary 
conditions in R . 

When n = 0, but m ^ 0, we set 



D m o(t)= / E(t,x + (2ma,0),x)dx 

JRo 
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— / / E(t,xi + 2ma,—X2,xi,X2)dxidx2 

J -a J 4\x 2 \ 2 <t 2 -4m 2 a 2 

— / / E(t,xi + 2ma,b — x 2 ,xi,b + x 2 )dxidx 2 - 

J -a-h\x 2 \ 2 <t 2 -4m 2 a 2 

As t [ t c , the domain of integration in the second and third integrals in D m0 (t) 
shrinks to {\xi\ < a, x 2 = 0}. Since these limiting cases correspond to the in- 
tegration of E(t,x,y) over {x = y + (2ma, 0), \yi\ < a, y 2 = 0} and {x = 
y + (2ma, 0), |yi| < a, y 2 = b} respectively, we will associate the second inte- 
gral with the side {0 < x\ < a, x 2 = 0} of the original rectangle R, and the third 
with the side {0 < xi < a, x 2 = b}. The first integral is paired with the interior of 
R. 

When m = 0, but n ^ 0, we set 



D 0n (t)= E(t,x + (2nb,0),x)dx 
JRo 

/ E(t, — xi, x 2 + 2nb, x\, x 2 )dx\dx 2 

J4x 2 <t 2 -4n 2 b 2 

E(t, a — X\, x 2 + 2nb, a + x\,x 2 )dx 1 dx 2 , 



-b J 4xl<t 2 -4n 2 b 2 

and we pair the first in integral with the interior of R, the second integral with the 
side {0 < x 2 < b, x 2 = 0} and the third with the side {0 < x 2 < 6, xi = a} 

Finally, when m = n = the term D 00 consists of the remaining terms in 
Tr(D(t)): the integral j Ra E(t,x,x)dx from e ++ paired with the interior R, four 
integrals from e" 1 and e ^ with m = and n = respectively, paired with 
the sides of R, and the four integrals from e in (11) which we pair with vertices 
(0,0), (a, 0), (0, b) and (a, b) of R. In this way we set up a one-to-one correspondence 
between the terms in D mn {t) and the k-simplices in R and dR corresponding to 
the zero components of (m, n), with k = 2, k = 2, 1 or k = 2, 1, depending on the 
number of zero components. This correspondence generalizes to higher dimensions 
(see Section 5). 

§3. Consequences of analyticity 

In this section we draw conclusions from the analyticity of the terms in the 
sums e ++ {t), e" 1 (i), e *~(t) and e (i). In [ER1] we showed that the terms 
J R E(t,d + x,x)dx from e ++ (i) were real analytic for \t\ > \d\ by appealing to 
results on analytic wave front sets. One can also prove the results we need on the 
other terms in the trace using analytic wave front sets. However, in this paper we 
will also use a more classical approach due to Hadamard [H] in the form presented 
by Hormander in [Ho]. 

Let E + (t,x,y) be the forward fundamental solution for the wave equation, i.e. 
the solution of ( - A + Q(x))E + = 5 {0 , y) in R t x K™, such that E + = for t < 0. 
Then E + can be written as 

oo 

E+(t,x,y) = ^a u (x,y)e u {t, \x - y\), (14) 

u=0 
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where 

e„(t, \x - y\) = 2 -2-i vr (i-")/2^+(i-»)/2^2 _ |x _ y|2) 

for t > 0, e„ = for t < 0. For a > — 1 the distribution is defined by = 
(r(a+l)) _1 s a for s > and <Y"(s) = for s < 0. One defines it for all a by analytic 
continuation in a; in particular, X+ k = <5( fc_1 ) for k £ N. The normalization of 
e„(t, | a; — y|) is chosen so that (d 2 — A)e l/ = ve„-i for v > , and eo(t, |x — y|) is 
the forward fundamental solution when Q = 0. The coefficients are determined 
by the recursion 

va v + (x - y) ■ d x a v + Qa v _ x - A x a v _ x = (15) 

or solving (15) 

a u (x,y) = / s u ~ 1 [A x a u - 1 (y + s(x-y),y)-Q(y + s(x-y))a u - 1 (y + s(x-y),y)}ds. 
Jo 

Starting with a (x, y) = 1 and using (15) repeatedly, one completes the construction 
of E + . The coefficients a\ and 02 are given by 



ai(x,y) = - Q(y + s(x - y))ds 
Jo 



(16) 



and 



02 



(x,y) = - [ s (l-s)AQ(y + s(x-y))ds + l([ Q(y + s(x - y))ds) 2 (17) 
■/ 1 Jo 



For general Q the sum in (14) is a singularity expansion: if one truncates the 
sum, then the difference of the fundamental solution and the truncation is more 
regular than the last term in the truncation. However, when Q is real-analytic 
one sees immediately that the a„'s are real-analytic. Moreover, by a majorization 
argument one can show that the sum in (14) is convergent for t sufficiently small as 
a power series in the variable (t 2 — \x— y\ 2 )+- Thus we have the following convergent 
expansions for E + near (t 2 — \x — y\ 2 ) + = 0. For n even the function 

A^(t x v) = V ^(x,y)(t 2 -\x-y\ 2 r 
A+[z,x, y) ^ 22 , +%(n _i)/ 2r(zy + (3 _ n)/2) 

is real-analytic, and we can write 

E + (t, x, y) = A + (t, x, y){t 2 - \x - y] 2 )^ 2 , (18) 

where (t 2 — \x — y\ 2 )+ n ^ 2 is interpreted in distribution sense as above. When n is 
odd, one keeps the first (n — l)/2 terms in (14) and combines the remaining terms 
as in (18). Let E_(t,x,y) = —E + (—t,x,y) and let E (t,x,y) be the odd extension 
of E + (t, x, y), i.e. E (t,x,y) = E + (t,x,y) — E_(t,x,y). Note that E (t,x,y) is 
the solution of (J^ — A + Q)u = in 1 x W 1 satisfying the initial conditions 
E o (0,x,y) = 0, ^£0(0, x,y) = S(x - y). Therefore E(t,x,y) = j^E (t,x,y) is the 
fundamental solution for the initial value problem in (2). 
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Using (18) one can give convergent expansions of the trace components D mn (t) 
near their singularities. Consider the expression (11). Assuming t > and making 
change of variables x = ty in each integral in (11), we obtain 

e ""(*) = ^ / , [A + (t,-ty,ty)+A + (t,(0,b)-ty,(0,b)+ty)+ 

J\v\<i/2 v 1 - 4 lyl 

+A+(t, (a,0)-ty, (a,0)+ty)+A + (t,(a,b)-ty, (a, 6) + ty)]dy) 

for < t < e. This gives as power series expansion for e (i) for i > 0. 

The same reasoning with the substitution x\ = (t 2 — 4n 2 6 2 ) 1 / 2 yi yields an ex- 
pansion for the n-th term in the summation for e '"(t) in (12) for t > \2nb\, and 
with the substitution for = (t 2 — 4 : n 2 a 2 ) 1 / 2 y 2 it yields an expansion for the n-th 
term in the summation for e H (t) in (13) for t > \2na\. Since solutions of the 
wave equation propagate at unit speed, these terms are zero for < t < \2nb\ and 
< t < \2na\ respectively. Finally, using (18) without substitution, one gets an 
expansion for j Rg E(t, d + x, x)dx is real-analytic for t > \d\. 

Using d = (2ma, 2nb), these arguments show that D mn (t) is equal to zero when 
\t\ < |d|,and that it is real-analytic for \d\ < \t\ < \d\ +e. However, we need to know 
that D mn (t) is real analytic on the half lines where \t\ > \d\. For this we will use 
well-known results on analytic wave front sets. In particular, the following wave 
front calculations are based on Theorems 3.10, 5.1 and 7.1 of [H62]. 

Using the coordinates (t, x, y, r, £, rj) on the cotangent bundle, the analytic wave 
front set of E satisfies 

WF A (E) C{x = y + ti/\Hr 2 = \C\ 2 ,V = -£}■ 

Since the normal bundle to {(t, x + d, x)} does not intersect WFa{E), E(t,x,x + d) 
is well-defined as a distribution and WFA(E(t,x,x + d)) C {t 2 = \d\ 2 }. Thus the 
contribution to D mn {t) from e ++ (i) is real analytic for \t\ ^ \d\. 

For the terms coming from e" 1 (t) we have to consider restrictions of E(t, x, y) 
to {(t, —x\,X2 + 2nb, xi, X2)} and {(t, a — xi, X2 + 2nb, a + xi, X2)}. Here again 
the normal bundles do not intersect WFa(E). Thus the restrictions E(t, —x±,X2 + 
2nb, xi, x 2 ) and E(t, a — x±, x 2 + 2nb, a + x\ , X2) are well-defined with analytic wave 
front sets contained in 

{(r,^,^) = (t,x,r / ,-e i + v'i^2+V2),r',e,v') €WF A (E)\ (t>x , y) }. 

Since the normal bundle to a line parallel to the xi-axis does not intersect this set 
when t 2 7^ 4n 2 6 2 , we conclude that the analytic wave front sets of 

/ E(t, — s, x 2 , s, x 2 + 2nb)ds and / E(t, a — s, X2, a + s, x 2 + 2nb)ds 
Jr Jr 

are contained in t 2 = 4n 2 b 2 . Hence, for n / 0, D 0n (t) is analytic for \t\ ^ \d\. 
Analogous arguments apply to the contributions from e~ + {t) and e (i), and this 
leads to the desired results for D m0 (t) and D 00 (t). 

Now we can use the analyticity results from the preceding paragraphs as in 
[ERT1] to find additional spectral invariants. We use the irrationality of a 2 /b 2 to 
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conclude that \d'\ = \d\ for d,d' € L implies that d' = (±e£i, ±d 2 )- We also assume 
for definiteness that a < b. Now we can recover individual terms D mn (t) in Tr D(t) 
as follows. 

All terms in Tr D(t) except D 00 (t) are zero for < t < 2a by domain of 
dependence. So from D(t) with t near zero we recover D 00 , and Tr D(t) minus 
Doo(t) is known for t > by analyticity. Continuing to reason in this way, new 
terms appear in the expansion of Tr D(t) at t = \d\ for each value of \d\. Since 
these terms are analytic for t > they are determined by the spectrum for all 
t> \d\. When both components of d are nonzero, the only terms come from e ++ (i), 
and we conclude that for d\d2 ^ 



^ / E(t, ±di + xi,±d 2 + X2,x)da 

4- J Ro 



± JRo 

is a spectral invariant. However, the symmetry Q{—X\,X2) = Q{x\,x 2 ) implies 
E(t,-x 1 ,x 2 ,-yi,y 2 ) = E(t,x,y). Thus 

/ E(t,-d 1 + xi,±d 2 + x 2 ,x)dx = J E(t,d 1 - x x ,±d 2 + x 2 ,-x 1 ,x 2 )dx, 
JRo JRo 

and making the change of variables X\ — > —X\, we have 

/ E(t,—di + x\,±d2 + X2,x)dx= / E(t,d\ + Xi,±d2 + X2,x)dx. 
J Ro J Ro 

Since we also have the symmetry Q(xi,—x 2 ) = Q(x\,X2), the same argument 
applies in the variable X2- Thus the contribution from e ++ (t) simplifies to 



4D mn (t) = 4 / E(t,d + x,x)dx. 

JRo 



(19) 



When d = (d\, 0), di 7^ 0, or d = (0, c^), ^2 7^ 0, we get additional contributions 
from e H and e h respectively. This gives the spectral invariants 

p poo pa 

D m o(t) = / E(t,d\ + x 1 ,x 2 ,x)dx 1 dx 2 - / (/ di + xi, -x 2 , x)dxi)dx 2 

</ -Ro J —00 J —a 

/oo /-a 

(/ di +xi,6-x 2 ,xi,6 + £ 2 )(ixi)<ix2 ( 20 ) 
-00 1/ — a 

for di = 2ma 7^ 0, and 

/» pb poo 

Don(t) = / E(t,x 1 ,d 2 + x 2 ,x)dx 1 dx 2 - / (/ -Xi, d 2 + ^2, x)dx 1 )dx 2 

J R J —b J — oo 

cb poo 



/o poo 
(/ a — xi,d2 +X2,a + x,X2)dxi)dx2 (21) 

./ — 00 

for d 2 = 2n6 7^ 0. In (20) and (21) we again made use of the observation that 
changing di to —di does not change D mn (t). 
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§4. Heat Traces 

Let D mn (t) be the wave trace spectral invariant from (19). Then 

—j=\ e- s2/4t D mn (s)ds = G(t,d + x,x)dx, 

where G(t, x, y) is the fundamental solution for the initial value problem 

u t = An - Qu in R + x R 2 , «(0, x) = f(x). 

Theorem 5.1 in [ERT1] gives the asymptotic behavior of G(t,x + Nd + e,y) for 
d, e G L as N — > oo. To describe this we need the "reduced potentials" q<i{x) = 
Jo Q( x + sd)ds and the associated fundamental solutions G d (t,x,y) for the initial 
value problems 

u t =Au- q d u in R + x R 2 , it(0, x) = f(x). 

Then 

m, x, ») - g„ (( , ,,,)!< | ;e -i.-»,"/«_i_ (1 + (22) 

where x' — y' denotes the component of x — y perpendicular to d and C is uniform in 
x and y on bounded intervals t 6 [0,T]. In [ERT1] the estimate (22) is only stated 
for specific choices of x and y, but the estimate for general x and y is proven. 

The structure of G^ makes (22) useful. Since Q(x + e) = Q(x) for e € L, 
inherits this property, and in addition q d is constant in the direction d, qd(x + sd) = 
qa{x) for s£R. It follows that Gd is the product of a free heat kernel and a heat 
kernel determined by qd, i.e. 

x = g^zM ^,,, rf , (23) 

where Hd(t, x + sd, y + rrf) = H d (t, x, y). Letting z denote signed distance from the 
origin on the line d ■ x = 0, Hd(t,z,w) is the fundamental solution for the initial 
value problem for u t = u zz — qd{z)u. Letting y = Nd + e + x in (22) and letting iV 
go to infinity, it follows that for e, d € L with tii^ 7^ J Rn G(t, x + iVd + e, x)dx, 
N € N, determines J R G^(x + e,x)dx for any d G L\0 and e G L. Hence the 
invariants in (19) determine 

/ G d (t,e + x,x)dx (24) 

■/Ho 

for each eel and d & L with di(f 2 7^ 0, and the Dirichlet spectrum for —A + q(x) 
on determines 

V / G d (t,e + x,x)dx. (25) 

This sum is the heat trace, ^2exp(—X n t), for —A + q d (x) on R with the periodic 
boundary condition, u(x + e) = u(x), e G L. Thus the Dirichlet spectrum for 
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—A + q(x) on R determines the periodic spectrum for —A + q d (x) on R for any 
d G L with d\d 2 7^ 0. 

If we choose 5 G L* such that 5 ■ d = and {<5 • e, e G L} = Z, then we can write 
the Fourier expansion of q d in the form 

00 00 
Qd(x) = ^ a ke 27Tik5 ' x , i.e. q d (x) = v(6 ■ x) where v(s) = ^ a k e 27Tiks . 

k = — oo k= — oo 

In §3(a) of [ERT] we showed that the periodic spectrum for —A + q d (x) on R 
determines the periodic spectrum for —\5\ 2 d 2 /ds 2 + v(s) on [0, 1] - and vice versa. 
Note further that Q(x) = Q(—x) and Q(x + d) = Q(x), d G L, imply that c/d(x) = 
Qd(—x). Hence = a_fc and = v(—s). The set of even potentials on [0,1] 
with the same periodic spectrum is either finite or a Cantor set (see [GT]). Thus, 
if —A + q and —A + q have the same Dirichlet spectrum on R, for did% 7^ the 
reduced potential qd must belong to the Cantor set determined by qd- 

The derivation of (25) applies equally well for d = (di,0), d\ 7^ 0, provided 
that e2 7^ 0, but in the remaining case, = 0, we need to use the invariants from 
(20). The same remark applies to d = (0,^2), d2 7^ with the invariants from (20) 
replaced by those in (21). We will only give the argument for d = (di, 0). The case 
d = (0,d 2 ) is completely symmetric. 

When d = (d\,0), q d {x) = ^2(^2) and Gd is given by 

G d (t,x,y) = -^=e-^-y^/ it H(t,x 2 ,y 2 ), 

where H is the fundmental solution for the initial value problem 

Ut = u X2X2 - q 2 (x 2 )u, n(0, x 2 ) = f(x 2 ). 

As in (19) the Dirichlet spectrum for —A + q on R determines f R G d (t, (0, e 2 ) + 

x,x)dx for e 2 7^ 0, and hence it determines $ h _ b H(t,x 2 + e 2 ,x 2 )dx 2 for e 2 7^ 0. 
Applying the asymptotic argument that lead to (25) to (20) with e = 0, we conclude 
that the Dirichlet spectrum on R also determines 

p pa poo 

I G d (t,x,x)dx- (/ G d (t,x 1 ,-x 2 ,x 1 ,x 2 )dx 2 )dx 1 

J Ro J —a J —00 

/a poo 
(/ G d (t,xi,b- x 2 ,xt,b + x 2 )dx 2 )dx 1 , 
-a J —00 

and hence it determines 

/b poo poo 

H(t,x 2 ,x 2 )dx 2 - / H(t, -x 2 ,x 2 )dx 2 - / H(t,b - x 2 ,b + x 2 )dx 2 . 
-b J — 00 J — 00 

Combining these two results, the Dirichlet spectrum on R determines 

00 „t 

y~] J H(t,2nb + x 2 ,x 2 )dx 2 ] 



n= — oo 
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/OO />00 
H(t, -x 2 ,x 2 )dx 2 - / H(t,b - x 2 ,b + x 2 )dx 2 . (26) 
-OO J — oo 

The arguments of Section 2 applied to the Dirichlet problem on [0, b] show that (26) 
is the heat trace for the Dirichlet problem for — (d/dx 2 ) 2 + q 2 (x 2 ) on [0, b]. In other 
words the Dirichlet spectrum for — A + q on R determines the Dirichlet spectrum for 
— (d/dx 2 ) 2 + q 2 (x 2 ) on [0,6] when d = (di,0). Analogously, the Dirichlet spectrum 
for — A + q on R determines the Dirichlet spectrum for — (d/dxi) 2 + qi(xi) on [0, a] 
when d = (0, d 2 ) 

§5. More General Boundary Conditions and Higher Dimensions 

The results of the preceding sections generalize to regions R = {x G K n : < 
Xi < (Xj, i = 0, ..,n} with either Dirichlet or Neumann boundary conditions on 
each face of dR. The arguments used earlier extend to this case, because (a) the 
"Method of Reflection" can be used to construct fundamental solutions, and (b) 
essential features of these fundamental solutions remain unchanged in the more 
general setting. We explain these points in what follows. 

In one space dimension one can check directly that the fundamental solutions 
for uu = u xx — q(x)u, u(0,x) = f(x), ut(0,x) = 0, on the interval [0, a] with 
Dirichlet or Neumann conditions on x = and x = a can be expressed in terms 
of the fundamendal solution, E(t,x,y), for the initial value problem u tt = u xx — 
Q(x)u, u(0,x) = f(x), u t (0,x) = 0, on the line, where Q is the extension of q to 
R as before, Q(x) = q{—x) for — a < x < and Q(x + 2na) = Q(x). The formulas 
are as follows: 

for the boundary conditions u(t, 0) = u(t, a) = the fundamental solution is 

oo 

E DD (t,x,y) = [T DD E](t,x,y) = ^ [E(t, 2ma + x, y) - E(t, 2ma - x, y)\ (27) 

m= — oo 

for the boundary conditions u(t, 0) = u x (t, a) = the fundamental solution is 

oo 

E DN (t,x,y) = [T DN E](t,x,y) = £ (-l) m [E(t, 2ma + x, y) - E(t, 2ma - x, y)] 

m= — oo 

(28) 

for the boundary conditions u x (t, 0) = u(t, a) = the fundamental solution is 

oo 

E ND {t,x,y) = [T ND E](t,x,y) = £ (-l) m [E(t, 2ma + x, y) + E(t, 2ma - x, y)} 

m= — oo 

(29) 

for the boundary conditions u x (t,0) = u x (t,a) = the fundamental solution is 

oo 

E NN (t,x,y) = [T NN E](t,x,y)= £ [E(t, 2ma + x, y) + E(t, 2ma - x, y)] (30) 

m= — oo 

These formulas can be checked as in §2. Note, for example, that in (23) we take 
an odd extension from [0, a] to [—a, a] to satisfy the Dirichlet boundary condition 
at x = and then we take an antiperiodic extension from [—a, a] to (— oo,+oo) to 
satisfy the Neumann boundary condition at x = a. 
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The linear operators T a p make it possible to express the fundamental solutions 
for boundary value problems in the region R above compactly. If we let denote 
T a p acting on the variable Xj, then the fundamental solution for u tt = Au — qu, 
u(0,x) = f{x), u t (0,x) = in R with Dirichlet or Neumann conditions on i ; = 
and Xi = ai, i = 1, n, is given by 

E afj (t,x,y) = [T^ i Tl P2 ---T^ n E](t,x,y), (31) 

where for each i, ai and Pi are either D or N. In (31) E(t, x, y) is the fundamental 
solution to the initial value problem u tt = Au — Qu, u(0,x) = f{x), ut(0,x) = 
in R n with Q(±xi, ±X2, ±x n ) = q(x) for x G R and Q{x\ + 2m,\a\, x n + 
2m n a n ) = Q(x). In the notation E a p in (31) we think of a and (5 as vectors. 
Formula (31) reduces to (2) when n = 2 and one imposes the Dirichlet condition 
on the whole boundary. 

iFvom (31) one sees that the only terms in the expansion of E a p(t,x,y) which 
do not have arguments of 2midi — Xi for some % are 

{-l) m ^ + - +m ^E{t,2m 1 a 1 +x 1 ,...,2m n a n + x n ,y), (32) 

(mi,..,m n )eZ™ 

where . . ,i r } is the set of indices i such that ^ fa. This is the fundamental 
solution for u u = Au — Qu in the domain R = {x : —ai < X{ < a^} with 
anti-periodic boundary conditions in the variables Xi with i = i\, ..,i r and periodic 
boundary conditions in the remaining variables. When we pass to traces as in 
(4), the contribution from these terms will correspond to the wave trace for this 
problem. The key observation is that all other terms in the traces will "telescope" 
as in §2 to integrals over R in the variables appearing in the form 277?,^ — X{. To 
see this note that the contributions of these terms to the trace when = /3j will 
be summations in the index m, of the form 

E(t, . . . , 2rriiai — Xi, . . . , x)dxi ■ ■ ■ , (33) 

rrii = — oo 

where the dots before and after the formula are intended to indicate summations 
on other indices and integrations in other variables. The reasoning which lead to 
(10) reduces (33) to 

/oo 
[E(t, ... , -Xi, ... ,x) + E(t, . . . ,ai - Xi, . . . ,ai + Xi, . . . )]dxi ■■■ . (34) 
-oo 

When «j t^ 1 Pi, formula (33) is replaced by 

oo „ Qi 

••• / E(t,...,2m i a i -x i ,...,x)dx i --- . (35) 

m i = — oo J —ai 

To see the effect of the factor (— l) m * recall the formula from Section 2 (rephrased 
in the current notation) 




E(t, 2midi — Xi, x)dxi = 
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/CO pOO 
E(t,...,2m i a i -x i ,...,x)dxi- / E(t, 2{rrn - 2)a, 
-ai J —ai 



Xi , . . . , X^dXi . 



The point here is that the cancellations arise separately between the terms corre- 
sponding to even indices m; and those corresponding to odd indices. Hence, the 
terms in (35) reduce to 



/: 



[E(t, -Xi, . . . , x) - E(t, ...,a,i-Xi,...,a,i+Xi,... )]dxi ■■■ . (36) 



Now we can apply the arguments from Section 3 and Section 4 in this more 
general setting. Very little modification is needed. We make the assumptions: 

A) the extension Q of q to M n described above is real-analytic on M n , and 

B) the numbers a\ , a\, are linearly independent over the rationals. 

We also introduce the rectangular lattice 

L = {(2m 1 a 1 ,2m 2 a 2 , ■■■,2m n a n ) : (mi,...,m n ) € Z n }, 

and its fundamental domain 

Ro = {x £ M. n : —ai < Xi < ai, i = 1, n}. 

Note that it follows from B) that \d\ = \d'\, d, d' € L is equivalent to = ±c^, i = 
l,...,n. 

As was done at the end of Section 2, the distribution wave trace 



2" / E af3 (t,x,x)dx= f E aP (t,x,x)dx= [ [T^T^ 2 ---T^ n E](t,x, 

J R J Re, <J Rn 



x)dx 



' Ro •* Ro 

can be represented as a sum 



I E a p(t,x,x)dx = J2 E U( t )> ( 37 ) 

where the term E^Jt) corresponds to d G L. As in Section 3 we see that each 
term E^Jt) is real analytic when \t\ > \d\ and equal to zero when \t\ < \d\. There- 
fore J2\d\=\d\ Eip(t) is a spectral inveriant for any d G L. When d has no zero 
components, 

E*p{t)= I E(t,d + x,x)dx. (38) 
■JRo 

When di l = • • • = di m = and d has no other zero components, the term E^ ^(t) 
is the sum of an integral of form (38), 2m integrals associated with the 2m faces 
x i± =0, x i2 = 0, ■ ■ ■ , x im =0, Xi t = a ix , ■ ■ ■ , x im = a im and integrals associated 
with (n — fe)-dimensional edges, 1 < k < m . For the uniformity of notation we 
shall consider a face as an (n — l)-dimensional edge and each vertex of R as a 
0-dimensional edge. Note that an (n — /c)-dimensional edge is the intersection of k 
faces. As in Section 2, all integrals except those of the form (38) tend to integrals 
over their associated edges when \t\ j \d\. 
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When d has no zero components, there are the 2 n terms in (32) corresponding 
to d' with \d'\ = \d\. Since the argument leading to (19) does not depend on the 
number of variables, each of these terms (up to a common factor of ±1 determined 
by the boundary conditions) equals (38). Thus, the wave trace (37) determines the 
partial traces (38) for each d 6 L with no zero components. 

Passing to heat traces as in Section 4, from (38) we conclude that, when Nd + e 
has no zero components, the wave trace determines 



/ G(t,Nd + e + x,x)dx. 



'Ro 

Taking the limit as N — > oo, we obtain 

/ Gd(t, x + e,x)dx, (39) 
JRo 

where Gd is the heat kernel associated with the potential qd given by 

q d (x) = / Q(x + sd)ds. 
Jo 

When it is possible to use this argument for all e 6 L, we conclude that the wave 
trace (37) determines 

V] / G d (t,e + x,x)dx. 

eeL jR " 

This is the trace of the fundamental solution for the heat equation on Rq with 
periodic boundary conditions, i.e. u(t, x + d) = u(t, x), d £ L, for the potential qd- 
At this point results of [ERT1] can be applied. Using the notation from Theorem 
1.1, the argument from pp. 668-9 of [ERT1] shows that for each 5 £ S such that 
S ■ d = the invariants in (39) determine the periodic spectrum of —A + Q§ (<5 ■ x) . 
As before the periodic spectrum of —A + Qs{5 ■ x) determines and is determined 
by the periodic spectrum of 

-\S\ 2 ^ + Qs(s) 

on [-1/2,1/2], [ERT1, Section 3(a)]. Also, Q s is even (Qs(-s) = Qs(s)), and the 
set of even potentials with a given periodic spectrum on [0, 1] is either finite or a 
Cantor set, [GT]. If S G S has at least two nonzero components, we can choose 
d 6 L with no zero components so that S ■ d = 0. Hence for any e 6 L the vector 
Nd + e will have no zero components for iV sufficiently large, and we can recover 
the invariants in (39) for all e £ L. Thus for each 5 G S with at least two nonzero 
components, the spectrum of — A + q on R with the given boundary conditions 
determines the periodic spectrum of — \5\ 2 d 2 /ds 2 + Qs{s) on [0, 1]. 

We are now left with the 5 € S with only one nonzero component, i.e. 5 = 
{2ai)~ l ei for some i, where ei,...,e n is the standard basis for W 1 . In this case 
S ■ d = implies di = 0, and we can only recover the terms in (39) for e 6 L with 
ei 7^ 0. To obtain the analog of the results in §4 in this case we need to consider 
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the sum of terms in the wave trace (37) with t c = \d\ when the i-th component of 
d is zero, but the other components are nonzero. This sum is given by 



(~l) p [/ E(t,d + x,x)dx + 
\d\=\d\ Jr ° 

+ei / [E(t,d+x—2xiei,x)+€2E(t,d+x+(ai-2xi)ei,x+aiei)]dx]M 

J{\xj\<CLj, j^i, \xi\<oo} 1 

where p, e\ and €2 are determined by the boundary conditions and are the same for 
all d. In fact t\ = 1 if CKj = N and ei = — 1 if oti = D. Likewise €2 = 1 when = 
and e2 = — 1 when a, 7^ Once again the terms in this sum are independent of 
the sign of the components of d, and thus each term in the sum is the same. Passing 
to the heat trace, applying the preceding with d = Nd + e where di = = and 
taking the limit as TV — > 00, we conclude that the wave trace in (32) determines 

/ G d (t,x + e,x)dx + (40) 
JRo 

+e 1 / [G d (t,e+x-2x i ei,x)+e 2 G d (t,e+x+(a i -2x i )ei,x+aie i )}dxM 

J{\x j \<a j , jjti, |a;i|<oo} § 

Now can apply a variation of the argument (from [ERT1, pp. 668-9]) used earlier. 
If n > 2, we choose a d' £ L with d\ = but no other nonzero components such 
that d and d' are linearly independent. Then, taking e = Nd' + e 1 with e\ 7^ in 
(39) and letting JV^oowe recover 

/ G dd '(t,x + e',x)dx, (41) 

where G dd i corresponds to the potential (q d ) d >- However, for e' with e\ = we 
make the same substitution in (40) and pass to the limit, getting 

/ G M >(t,x + e',x)dx+ (42) 
JRo 



ei 



/ [Gdd' (t, e'+x-2xiei,x)+e 2 G dd >(t, e / +x+(a i -2x i )e i , x+a^ei)]^.! 

J{\xj\<aj, jjii, |xi|<oo} U 



We continue this until (((q d ) d ')d")--- = Qs(^ ■ x) for d = a i 1 e i . Letting G5 denote 
the heat kernel associated with the potential Q§(5 • x) = Q$(xi/ai), we have 

G s (t,x,y) = (47rt)-(™- 1 )/ 2 e -l-'-f'l 2 / 4t ft (t,x l ,y l ) 

where x' and y' have the i-th components omitted, and gi is the fundamental 
solution for the initial value problem for u t = d 2 u/dx^ — Qs(xi/ai)u. Thus the 
invariants in (36) and (37) are equivalent to 



/ai 
9i(t,Xi + ei,Xi)dxi and 
-a,; 



19 




/ gi(t,Xi,Xi)dxi + ei / [gi(t, -Xi,Xi) + e 2 ft(i, a; - Xi, + 



J —a t J -co 



Thus, as in the argument at the end of Section 4, undoing the simplification from 
the telescoping terms, these invariants determine 



Since (43) is the heat trace corresponding to the boundary condition on 
[0, ai], we have the generalization of the result of Section 4: 

Theorem 5.1. Suppose that q and q are isospectral potentials on R such that 
the corresponding potentials Q and Q on M n are real-analytic, and suppose that 
{a 2 ,--- , a^} are linearly independent over Q. Then for each 8 6 S with more 
than one nonzero component — \S\ 2 d 2 /ds 2 + Qs(s) and — \8\ 2 d 2 /ds 2 + Qs(s) have 
the same periodic spectrum on [—1/2,1/2]. For 5 = (2ai)~ 1 6i the operators 
—d 2 /dx 2 + Qs{xi/ai) and d 2 /dx 2 + Qs{xi/ai) have the same spectrum for the 
boundary condition (aj,/3j) on [0, <Zj], z = l,...,n. 

Remark 5.2. When Qs = for all 5 6 S with more than one nonzero component, 
Theorem 5.1 gives all the constraints imposed by isospectrality: when the operators 
—d 2 /dx 2 + Qs{xi/ai) and d 2 /dx 2 + Q&{xi/ai) have the same spectrum for the 
boundary condition (a^,/%) on [0, ctj], i = l,...,n, q and g are isospectral. This 
is not what one expects when other directional potentials are nonzero. In [ERT2] 
we used higher order terms in the asymptotics of G(t, Nd+e+x, y) to get additional 
contraints on isospectral potential and establish rigidity modulo lattice isometries 
in some cases. We plan to carry this out for the boundary conditions considered 
here in the future. 

As noted in the Introduction, the set of Qs on [—1/2, 1/2] with a given periodic 
spectrum is necessarily discrete since the Q^'s are even. When one works with 
general L 2 potentials, this is not true for the operators —d 2 jdx 2 + Qs(xi/ai) on 
[0, a^]. See [PT, Chpt. 6] for (explicit!) examples of infinite dimensional manifolds 
of potentials with a given Dirichlet spectrum. However, when one restricts the 
admissible potentials to those considered here, the set of isospectral potentials for 
these operators is again discrete. 

Lemma 5.3. Assume that the potential Q on M is real-analytic and satisfies 
Q(x) = Q(—x) and Q(x + 2a) = Q(x). Then the spectrum of d 2 jdx 2 + Q(x) on 
[0, a] with either Dirichlet or Neumann boundary conditions at x = and x = a 
determines the periodic spectrum of d 2 /dx 2 + Q(x) on [—a, a]. 

Proof. The reasoning given earlier, specialized to n = 1, shows that the spectrum 
determines 




(43) 




and 




/ E(t, x,x)dx + ei / [E(t,-x,x) + e 2 E(t, 



a — x,a + x)]dx, 



(44) 



J — a J — oo 



where e\ equals 1 or -1 according to whether Neumann or Dirichlet conditions 
are imposed at x = 0, and e 2 equals 1 or -1 according to whether the boundary 
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conditions at x = and x = a are the same or different. Our goal is to show that 
as long as Q is a real-analytic potential satisfying Q(—x) = Q(x) and Q(x + 2a) = 
Q(x), the spectrum determines 

pa 

E(t,x,x)dx, (45) 

and hence determines the periodic spectrum. 

In one dimension the Hadamard formulas become 

E(t,x,y) = d t E (t,x,y), E (t,x,y) = E + (t,x,y) - E + (-t,x,y), and 

OO 

E+(t,x,y) = ^a v (x,y)e v (t, \x - y\), where e„(t,\x-y\) = 
1 (t 2 - \x- y\ 2 f + for t > 0, and E + (t,x,y) = for t < 0. 



2 2 "+ 1 r(i/ + 1 
Hence 



1 00 />a 

E(t, x, x)dx = ^ 2 2v+ir(v + 1) J a »(. x , x )dx[d t (t 2u sgn(t))]. 

We are going to show that the other terms in (44) contribute only even powers of 
t when one expands (44) in powers of t, and hence one can recover (45) from (44). 
We have (for t > 0) 

/OO 00 j />oo 

E(t,-x,x)dx = J^ —^—^ J aA-z,x)[d,{t 2 - 4x%]dx 

Note that in removing the term corresponding to v = in the last formula we used 
ao(x,y) = 1. Similarly, 

,00 oo^ 2 ^ ,1/2 

/ E(t,a-x,a+x)dx = V ^tTp? — T~n / a u (a+ty,a-ty)t 2 " '(l-4y 2 ) !y-1 dy 

7-oo „ =1 Z 1 1 17 + X J 7-1/2 

So to eliminate odd powers of t the Maclaurin series for a u (—x,x) and a v {a — 
x, a + x) must contain only even powers of x. From (16) 

ai (bx,x) = - I' Q((l + (b-l)s)x)ds = V- (1 + (& " 1) ^ Q( fc )(0)d 

which contains no odd powers because Q( fe )(0) = for k odd, and 

ai(a + bx, a + x) = — / Q(a + (1 + (b — l)s)x)ds 
Jo 



sx k 
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which also contains no odd powers because Q^ k \a) = for k odd. From the 
recursion relation (17) we have 




Jo 





Jo 



Using the induction hypothesis: a M (6x, x) has an expansion in even powers of x for 
fi < u, we can conclude that a u (bx,x) has an expansion in even powers of x (note 
that d^.a u -i(bx,x) = x~ 2 d 2 a„-i(bx,x)). Likewise 



a u (a + bx, a + x) = I s u 1 [d 2 a u -i(a + (1 + (6 — l)s)x, a + x) 
Jo 

-Q(a + (1 + (b - \)s)x)a v -i(a + (1 + (b - l)s)x, a + x)]ds. 

Hence, using the induction hypothesis: a fl (a + bx,a + x) has an expansion in even 
powers of x for ji < u, we can conclude that a„(a + bx,a + x) has an expansion in 
even powers of x. We complete the proof that the Maclaurin series for 
and a u (a — x, a + x) must contain only even powers of x by taking b = — 1 . Thus we 
have shown that the the spectral invariant in (44) determines the integral in (45) 
which combined with the spectral invariants preceding (44) determines the periodic 
spectrum. □ 

Lemma 5.2 combined with Theorem 5.1 gives the following generalization of 
Theorem 1.2. 

Theorem 5.4. If {a\ , ■ ■ ■ , a^} are linearly independent over Q, then for any choice 
a(3 of boundary conditions on the faces of R, any continuous curve q t of isospectral 
potentials such that the extensions Q t are real analytic for all t must be constant, 
Qt — Qo f° r an t- 

§6. Singularities of the wave trace and applications. 

In this section we drop the requirements that Q(x) is real-analytic and {a\ , a^} 
is linearly independent over Q, and assume only that Q(x) 6 C°°(lR n ). We shall 
the study of the singularities of Tr E a p(t). We assume that d € L\0 and find the 
singularities of Tr E a p at t = |<Z|. The case d = is discussed in [ER]. Using the 
notation of §5, the only terms of Tr E a/ s that are singular at t = \d\ are 



All other terms in Tr E a p are C°° near t = \d\. We shall find the singularity 
expansion of (46) at t = \d\, and the terms of this expansion will be spectral 
invariants. 





(46) 



\d\=\d\ 
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It follows from the Hadamard construction (see (14), (16), (17) ) that the sin- 
gularity expansion of 

E(t,x + d,x)dx (47) 



Ir, 



Ro 

at t = \d\ has the form 



70 f R idx^x7 L {t 2 -\d\ 2 ) + ~ h f R q {x)dx^xf + \t 2 -\d\ 2 ) 

+ 72 j R q l( x )dx^- t xf +2 {t 2 -\d\ 2 ) + -.., (48) 

where ji, < i < 2, are constants independent of R and Q(x), and qd{x) = 
j Q Q(x + sd)ds is the reduced potential. In deriving (48), we noted that Q(x) is 
periodic and smooth, and therefore J R Q(x + sd)dx = j R Q{x)dx = 2 n j R q(x)dx 
and 

/ AQ(x + sd)dx = 0. 
jRa 

When d 6 L has no zero component, E^^t) consists only of the integral (47) (up 
to factor ±1 depending on boundary conditions). Now suppose that d^ = ■ ■ ■ = 
di m = and d has no other zero components. Then, as in Section 5, E^g{t) in 
addition to (47) contains a sum of integrals associated with edges of dimensions 
n — k, 1 < k < m. Consider, for example, an integral I c '(t), associated with the 
edge r c / = {x : x ix = c ix , •, x ik = c ik }, where the j-th component of c' = (c ix , c ik ) 
is either or a^.. We write x = (x',x"). This, as in Section 5, is a slight abuse 
of notation: x' = (x^ , Xi 2 , ■ ■ ■ ,Xi k ) denotes a subset of components of x - not 
necessarily the first k components. Near t = \d\ = \d"\ (up to factor ±1) 

I c ,( t ) = ^-[ [ E + (t,c' -x',x" + d",c' + x',x")dx'dx". (49) 
Jr c , JR k 

Substituting the expansion in (14) into (49), we get the expansion 

oo 

/ 2 c'(t) = ]T/ c> (t), (50) 



v=0 



where 



d 



l-n 



/^(t)=7^ J J k X~^(t 2 -\d\ 2 -Mx'\ 2 )a u (c'-x', X "+d",c'+x',x")dx'dx" 

(51) 

Making the change of variables x' = t^y' where ta = (t 2 — \d\ 2 )+ as in Section 3, 
expanding a v (c' — t^y', x" + d" , c' + t^y' , x") in a Taylor series in tay', and noting 

that x+ 2 V (s) is homogenous of degree + v, we get 

h'At) ^J2 b ^PQ- t xT +U+i+P (.t 2 - \d\ 2 ), (52) 

p=0 
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where 

bc>vO = 7vo([ xT^il-^vVW) I a„(c',x" + d",c',x")dx" 

JR k JY C , 

and the b c i vp are integrals over F c > of derivatives of a v {c' — x' , x" + d" , c' + x' , x") 
in x' at x' = 0. Using (15) we see that 6 C /„ P is a sum of integrals of polynomials in 
Q{x) and its derivatives over T c i. Note that integrals in I c >,v{t) having odd powers 
of y' in the Taylor formula will be zero. Note also that formula (52) remains true 

1 — 71 | ^ 

when x+ 2 is a distribution, i.e. when + v < — 1. To check this we use the 
representation 

X? + V " \d\ 2 ~ 4|*f) = ^ t J t ) N xT + " +N {t 2 -W- 4|*f), 

where +u + N > — 1. Now we can change variables x' = (t 2 — \d\ 2 )+y' as above 
and then apply (^-^) N to get (52). Collecting (6), (50), and (52), we get 

E - E E E ^(^Ix?^* 4 *^ - \d\ 2 ) (53) 

|<2| = |d| |rf| = |rf| fc=1 ^P=° 

Therefore ^|J|=|d| S^+ P +^=r b v k P {d) are spectral invariants for r = 0, |, 1, |, • • • . 

Here we will briefly describe the simplest of these spectral invariants. Since the 
coefficient ao is identically 1, b v k p (d) with v = does not depend on q and it is 
zero for p > 0. Thus terms depending on q first appear for r = 1. The invariant 
for r = 1 is the sum of &ioo(^) an d &o2o(rf)> an d it has the form ci J R q(x)dx + c%, 
where c\ only depends on n and C2 is determined by R. Thus J R q{x)dx is a spectral 
invariant as claimed in Section 1. 

For r = 3/2 the terms involving g appear when d has one or more zero compo- 
nents and correspond to (u, k,p) = (1, 1, 0). Hence they correspond to integrals of 
q over the edges of R of codimension one. Other contributions to the invariant for 
r = 3/2 appear only when d has three or more zero components and correspond to 
(u,k,p) = (0,3,0). 

The invariants needed for the proof of Theorem 1.3 correspond to r = 2. When 
r = 2, contributions to the invariant come from (v,k,p) = (2,0,0), (1,2,0) and 
(0,4,0). The (2, 0, 0)-term is the third term in (48), but the (1, 2, 0)-term also 
depends on q, since it is a linear combination of the integrals of q over edges of 
codimension 2. Since the index k must be less than or equal to the number of zero 
components of d, we see that 

£ / {q d {x)fdx (54) 
\d\=\d\ jR 

with d 7^ will be a spectral invariant if n = 2, but it will not be an invariant for 
n > 2. 

Proof of Theorem 1.3. Suppose without loss of generality that a < b. We 
consider the case a < b first. Take v\ = (a, 0). Then d = (±a,0) if \d\ = \vi\. 
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Therefore, using the invariance of (54), f Ro q 2 i (x)dx = j Rg q 2 1 {x)dx. Take v 2 = 

(0,6). If \d\ = |f 2 1 then either d = (0, ±6) or d = kv\ if b = ka, k € Z. In both 
cases we have 

I q 2 V2 (x)dx + C f q 2 Vi (x)dx= f q 2 2 (x)dx + C f q 2 Vi dx, 

J Rq J Rq J Rq J Rq 

with C = in the first case and C = k 2 in the second. Therefore 

/ q% 2 (x)dx = / ql 2 (x)dx. 

J Rq J Rq 

In the case when a = b we have by the invariance of (54) 

/ {q 2 Vl + Q% 2 )dx = / {q 2 Vl + ql 2 (x))dx. 

JRo J Rq 

Take an arbitrary d € L, such that d\d% ^ 0, and consider all d G L such 
that \d\ = \d\. These d's may have the form kxvi or k 2 V2 or they may have both 
components nonzero. We label the ones with both components nonzero d' . We 
have by (54) 

k\ / q 2 Vi dx + kl / q 2 2 dx+ / q%( x ) dx 

JRo JRo \d'\=\d\ jRo 




Since qj,(x) = for all d', (55) implies q^i{x) = for all d'. This implies q = 
Qi{xi) + q 2 {x 2 ). □ 

In case when the lattice L satisfies the conditions of Theorem 5.1 we can easily 
prove a more general result. 

Theorem 6.1. Suppose q(x) and q(x) have the same a(3 spectrum on R and their 
extensions Q(x) and Q(x) are sufficiently smooth in R n . If {a 2 , • • • , a^} is linearly 
independent over Q, q(x) = <?i(xi)H Vq n {x n ) implies q(x) = qi(x\)-\ Vq n {x n ) 

Proof: Suppose d = (di,--- ,d n ) has no zero components. Since we assume 
{a 2 , • • • , a^} is linearly independent over Q, |J| = \d\ if and only if d = (±di, ±d 2 , • • • , ±d n ).( 
Since none of these d's have zero components, (47), summed over {d : |d| = \d\} 
is a spectral invariant. Thus from (48) we see that (54) is a spectral invariant. 
Since q(x) = qi(x) + ■ ■ ■ + q n (x) and d has no zero components, we have qa{x) = 0. 
Therefore by (54) qaix) = for all d with no zero components. If Qs were nonzero 
for some S with more than one nonzero component, we could choose a do with no 
nonzero components such that 5 ■ do = 0. Since this implies q~d 0> we have a 
contradiction. Hence Qs = whenever 8 has more than one nonzero component, 
and we have q(x) = qi(x\) + • • • + q n (x). □ 

Note that in the case of periodic spectrum the analog of Theorem 6.2 was proven 
in [ERT1] and [GK] under much weaker conditions on the lattice. 
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§7. Heat trace asymptotics when t — > 0. 

Consider the case of the Dirichlet boundary conditions in the rectangle R, n = 2. 
For \t\ small we have 

4 Tr D(t) = A)o(t) = 777 / E (t,x,x)dx 

/ / (£ (i,Xi, -X 2 ,Xi,X 2 ) + i?o(^^l,^-^2,^l,^ + ^2))d^l^2 

J-a Jix 2 2 <t 2 

- / (£ (t, -xi,x 2 ,xi,x 2 ) + £ (t,a - xi,x 2 , a + xi,x 2 ))dxidx 2 + e (i) , 

J-bUx\<V 

(56) 

where e (i) has form (11). Note that E (t,x,y) is an odd distribution in t: 



d_ 

di 



v=0 

where 



E (t,x,y) ^^2a u (x,y)e^\t,\x - y\), 
e i 0) (*, N) = e v(t, M) ~ e v(t, M), 



e v (t, \x\) = for t > 0, e u (t, \x\) = e u (—t, \x\) for t < 0. Note that e^\t, \x\) is a 
distribution in t depending smoothly on \x\. is odd in t and homogeneous of 
degree —1 + 2v when n = 2. In particular when \x\ = we have 



Note that P-v.\ = ^ ln|t| is the only odd distribution (up to a constant) homoge- 
neous of degree -1. 

Using (14), (16), (17), we recover (48) in the form 

f Id 2 If t 2 f 

/ E(t,x,x)dx = — \n\t\Aab / Q(x)dx + / Q 2 (x)dx + 0(t 4 ). 

JRo 27r dt 47r JRo 167r JRo 



Consider one of four integrals in (56) associated with sides of R, for example, 
§lh(t), where 

h(t) = / / E (t,a- x 1 ,x 2 ,a + xi,x 2 )dx 1 dx 2 . 

J-b Jixl<t 2 

Substituting (14), (16), (17) and making the change of variables x\ = ty\ we get 
for t > 

m _ f b f \ 1 1 t 2 (l-4y 2 y/ 2 f 1 



f f 1 1 t 2 (l — Ay 2 ) 1 ' 2 f 

/ / TTn a 2M/2 a — 1 / Q(a + tyi-2st yi ,x 2 )ds 



2 

Is 
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(1-4^)3/2 f (i- s ) s AQ(a + ty 1 -2sty 1 ,x 2 )ds 
Jo 



d yi dx 2 + 0(t 6 ). 



We expand Q(a + — 2s), x 2 ) in a Taylor series in iyi. Noting that integrals 
containing odd powers of y\ are equal to zero and extending I\ (t) for t < as an 
odd function we get 



h(t) = sgn t 



,6 

7i26 + t 2 7 2 y Q(a,x 2 )dx 2 



\ I ( 3 2 (°' X2 ) dx2 + 74 y Qx 2 1 ( a , x 2)dx 2 



+7 5 y AQ(a,x 2 )(fx 2 ) 



+ 0(t 6 ). 



Since Q(xi,x 2 ) is periodic in x 2 we have J_ b Q x 2 (a, x 2 )dx 2 = 0. 

Next we compute the asymptotic expansion of the four terms in e (t) associated 
with vertices of R. Consider, for example, Mrl 2 {t) where 



h{t) 



E (t, -xi,b - x 2 ,xi, b + x 2 )dxidx 2 . 



U{x\+xl)<t 

Making the change of variables x\ = ty\,x 2 = ty 2 , for t > we have 

tdy t 3 (l-4|y| 2 )5 f 1 



h{t) 



4|y| 2 <l 



27r(l-4|y| 2 )4 



4-7T 



Q(*2/i - 2isyi, 6 + ij/i - 2tsy 2 )ds 



dy+0{t l 



1 



= -f 6 t + 77 t 3 Q{0,b) + O{t 5 ). 
Combining the computations for all terms in (56), and evaluating the 7j's, we get 



/ D(t,x,x)d. 
Jr 



+ 



16 



L qds+t2 (wj/ {x)dx -kt^. 



-w 



2 7 



9/V 



(57) 



Here =06 is the area of R. \dR\ = 2a +26 is the length of the boundary dR, J^q 
is the normal derivative of q on OR, Pi, 1 < i < 4, are the vertices of -R. 

Denote by G(i, x, y) the heat kernel in R corresponding to Dirichlet boundary 
conditions. As in Section 4 we use the formula 



G(t,x,y) 



1 



2Virt 



D(s,x,y)e s ^ds. 



(58) 
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Applying this formula to (57) we have 



eft...)* = £w - ^|9bi + (i -ij R ^)ix) + SL 5 * 



9 2 



(59) 



gds )+0(t 2 



I 



To illustrate the flexibility of the method we are using we will do one more 
example. Consider the case of the Dirichlet boundary conditions when n = 3, R = 
{0 < Xi < cii, i = 1,2,3}. Denote by F c . ,i = 1,2,3, the face Xi = Ci of R where Cj 
is either or a*. Let T CiC . be the edge, T CiCj = F Ci nF Cj , i ^ j. The trace 2 3 D 00 (t) 
is a sum of the following terms: 

J Rg E(t, x, x)dx associated with the interior of R, 6 integrals associated with faces 

T Ci , 12 integrals associated with edges F ClC . and 8 integrals associated with vertices 
p 

Ci Cj Cfc 

Analogously to (57) we get 



f^D(t,x,x)dx = l0 5'\t)\R\+ ll j t p.v.- t Y J \^c l \+m 721 J q{x)dx + 722 I 1 ^ 



+ 



Ci ,Cj ,Cfc 



731 X / 1 ds+ ^32 ^ \ Pc i c i c k\ +1*1 741 / q 2 (x)dx + 742 5^ / qd 



/" f d 2 

751 XI / 1 2ds + 752 / o~2 ^ S + 753 



C i Cj Cfc j 



Ci -jCj , Cfc 



+ 0(i 3 ). (60) 



Here |r c J is the area of r ci ,|r ciCi | is the length of T CiCj and \P CiCj c k \ = 1 by 
definition. Using (58) we can easily obtain the heat trace expansion from (60). The 
derivation of (60) is the same as that of (57). We only mention a precaution needed 
to handle the change of variables in 



r C i jr 



e^\t, 2\xi\)dxidx2dx3 



since e^\t,2\xi\) is a distribution. We have 



e (t,2|xi|) = — 5(t 2 -4x 2 ) 
Since is odd in t, we have for (j) £ Cg°(R 1 ) 



(^°\t,2\ Xl \U) = (e (t,2\x\)A(t)-4>(-t)) = ^ (t8{t 2 -Ax\) m / ( t] 
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Note that 



4>(t)-4>(-t) 
t 



€ and t5{t 2 -Ax\) = ±£ t 6(t 2 -4xj) where 6(s) = 1 when 
s > 0,9(s) = when s < 0. Integrating 6(t 2 — 4xf) over x\ and changing variables 
X\ = tyi we get 



{ef\t,2\x 1 \),^)dx 1 



4tt 



t. 



dt t 



1 

47T 



0(t), 



Ht) - H-t) 



47T 



Similar arguments show that (e^\t, 0), <f(t)) = ^^'(0) = —^{S',<ft) 

With a little more effort one can compute the coefficients in (60). They are 
7o = -(27T)" 1 , 7l = -(Svr)- 1 , 721 = -(4vr)-\ 722 = (16)" 1 , 731 = (levr)" 1 , 
732 = -(64)" 1 , 741 = (327T)- 1 , 751 = -(64)" 1 , 752 = (1287T)" 1 and 753 = (128)" 1 . 
Applying (58) to the expansion in (60), we get the expansion of the heat trace 



-Mi* _ 



+ 



8(7ri) 3 / 2|i?l 16^^ |rCll + tV2 



8^// (x)dX+ 32^El r - 



1/2 



+ t 



^E/ r <^+^E/ r E 

Ci 1 c i Ci 1/1 =i Ci,Cj,C k 



+ o(t 3 / 2 ). 



These computations can, of course, be carried out for arbitrary assignments of 
Dirichlet and Neumann conditions on the sides T Ci . The only changes are as follows. 
The integrals in Dqq corresponding to T Ci with the Dirichlet condition are preceded 
by minus signs while those corresponding to F Cz are preceded by plus signs. The 
integrals corresponding to edges T c . c . with the same boundary condition on T Ct 
and Y c . are preceded by plus signs, and the others have minus sign. Finally the 
integrals corresponding to vertices are preceded by (— l) m where m is the number 
of adjoining faces with the Dirichlet condition. 
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